Astronomy & Astrophysics manuscript no. paperl4276 


©ESO 2010 


December 15, 2010 





Electron fraction constraints based on Nuclear Statistical 
Equilibrium with beta equilibrium 

A. Arcones 1 ' 2 , G. Martmez-Pinedo 2 , L. F. Roberts 3 , and S. E. Woosley 3 



1 Institut fur Kernphysik, Technische Universitat Darmstadt, Schlossgartenstr. 9, D-64289 Darmstadt, Germany 

2 GSI Helmholtzzentrum fur Schwerionenforschung GmbH, Planckstr. 1, D-64291 Darmstadt, Germany 

3 Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064 



o 

(N 
O 

Q 



in 
6 

> 

in 

oo 
m 

(N 
O 
O 



X 



received; accepted 



ABSTRACT 



The electron-to-nucleon ratio or electron fraction is a key parameter in many astrophysical studies. Its value is determined by weak- 
interaction rates that are based on theoretical calculations subject to several nuclear physics uncertainties. Consequently, it is important 
to have a model independent way of constraining the electron fraction value in different astrophysical environments. Here we show that 
nuclear statistical equilibrium combined with beta equilibrium can provide such a constraint. We test the validity of this approximation 
in presupernova models and give lower limits for the electron fraction in type la supernova and accretion-induced collapse. 
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1. Introduction 

Weak interactions are known to play a critical role in the struc- 
ture and evol ution of stars and the kind of nucleosy nthesis they 
produce (e.g. [Langanke & Martlnez-Pinedq[ l2003h . For nuclei 
near the valley of beta-stability, the relevant lifetimes of unsta- 
ble nuclei are well determined in the laboratory. However at the 
high temperatures and densities found in stars and stellar explo- 
sions, nuclei may be produced that are either very proton- or 
neutron-rich. Usually the rates for weak interactions involving 
such nuclei are only available from theory. Further, these nuclei 
exist in a distribution of excited levels, making even a theoretical 
calculation of their weak lifetimes problematic. 

Over the years, different groups have calculated weak- 
inter action rates for astrophysica l ap plications (for a review 
see lLanganke & Martinez-Pineda |2003l) . In the 1 980s Fuller, 
Fowler. & Newman dFuller et all 1 1980. Il982blal I1985L here- 
after FFN) calculated rates for electron capture, positron 
capture, beta decay, and positron emission for astrophysi- 
cally relevant nuclei. Their tabulations were based on an 
examination of all available information at that time and 
became the standard in the field. With the improvement 
of methods and computers, complete she ll-model calcula- 
tions were performed f or «/-shell nuclei dOda et all | l994) 
and for pf-shell nuclei (|La nganke & Martfnez-Pinedoi [2000; 
lLanganke & Martmez-Pinedoll200lL hereafter LMP). Although 
these new shell-model rates and the FFN rates agree rather well 
for sof-shell nuclei, there are appreciable difference at higher 
mass. The electr on-capture rates, and to a lesse r extent, the beta- 
decay rates (see Martm ez-Pinedo et all l2000h are smaller than 
the FFN rates by approximately one order of magnitude. 

In this paper we discuss an interesting limiting case that can 
be used to test existing rate tabulations and their implementation 
as well as to calculate astrophysical models in regions where 
rate information may currently be inadequate. Dynamic beta- 
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2009) occurs when sufficient time elapses at a given temperature 
and density for the electron mole number, Y e , to approach a con- 
stant. That is, for an ensemble of nuclei (but not necessarily indi- 
vidual nuclei), electron capture and positron emission occur at a 
rate balanced by electron emission and positron capture. While 
physically not the same situation as when neutrino absorption 
balances neutrino emission (thermal weak equilibrium), the so- 
lution to the abundance equations is the same in many of the 
astrophysical environments of interest, as we will show. Where 
this approximation holds, a set of nuclear abundances can be ob- 
tained, and thus a steady state value for Y e , that depends only on 
nuclear bulk properties (like binding energy and partition func- 
tio n), the temperature, and the den sity. Similar work was done 
bv llmshennik & Chechetkinl (1 197 lh . but our results are based on 
improved nuclear data present a broader overview of the possi- 
ble astrophysical scenarios where this approximation becomes 
useful. 

There are only a few places in nature where thermal weak 
equilibrium actually occurs, but they are interesting. Obviously 
the Big Bang is one, but there the abundances are simple, 
just neutrons a nd protons. Dynamic beta-equili brium exists in 
neutron star£] dBavm et all 1197 It IWeberi 1 1999b and may oc- 
cur briefly in the cores of massive stars after s ilicon burning 
( Aufderhe ide et all [l994bllat iHeger et all l2001allbl) . It may also 
happen in white dwarfs th at ignite carbon runaways at very high 
density dCanal et al.Lll990b . In the latter case, there is thought to 
exist a critical ignition density above which the reduction in elec- 
tron density by electron capture more than offsets the rise in ther- 
mal pressure due to the burning. Above this density whose exact 
value depends on composition, initial mass, and accretion rate 
dCanal et allll990h . it is expected that the white dwarf collapses 
to a neutron star rather than exploding as a supernova. This sce- 
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nario is still only theoretical since it was not yet observed. Just 
short of this critical density, which m ay be ~ 8 x 10 y g cirT 3 
dTimmes & Wooslevl 119921: IWooslevl 1 19971) . the star does ex- 
plode but ejects neutron-rich isoto pes (e.g., 48 Ca) that are diffi- 
cult to make elsewhere in nature (IWooslevl fl997). In addition, 
as the Chandrasekhar mass scales as Men K Y 2 , if a flame re- 
duces the Y e before the star begins to expand vigorously, it will 
collapse. The race between expansion and electron capture must 
be determined by hydrodynamical simulations with physical fi- 
delity, preferably in 2D or 3D. Until they are done, it is difficult 
to make quantitative statements, but it is expected that a 10% 
change in the terminal value of Y e would cause an approximately 
10% change in the terminal density and might affect whether 
e.g., the critical ignition den sity for collapse is 8 or 9 x 10 9 g/cm 3 
(Tim mes & Wooslevl 1 1992b . Instabilities could make the depen- 
dence non-linear. Thus it is important to know what the mini- 
mum achievable value is for Y e at a given temperature and den- 
sity. 

To proceed, we calculate the electron fraction and the com- 
position assuming thermal weak equilibrium for temperatures 
and densities that are relevant to various astrophysical environ- 
ments. Thermal weak equilibrium will only exist at such high 
temperatures and densities that nuclear statistical equilibrium 
(NSE) also prevails, and we assume that to be the case. We 
show that for temperatures, T < 10 GK, thermal weak equi- 
li brium reduces approximatel y to the dynamic beta-equilibrium 
of lTsuruta & Cameronl dl965). This allows to obtain results that 
depend only upon binding energies and partition functions, not 
individual weak interaction rates. At higher temperatures, ther- 
mal neutrinos are absorbed with a rate comparable to the elec- 
tron emission rate. When this happens thermal weak equilibrium 
and dynamic beta-equilibrium become different. Nevertheless, 
the equilibrium value of the electron fraction for dynamic beta- 
equilibrium can still be obtained using a rather simple descrip- 
tion of the relevant weak interaction rates. 

Our paper is organized as follows. We start with a brief 
review of NSE, thermal weak equilibrium and dynamic beta- 
equilibrium in Sect. [2] The results obtained from combining 
these assumptions are presented in Sect. [3] for broad ranges of 
temperature and density. In Sect. [4] we discuss the astrophysical 
implications of our results for presupernova models (Sect. 14. j} , 
proto-neutron star evolution (Sect. 14. 2I >. and accretion-induced 
collapse and type la supernovae (Sect. l4~3l l. Finally, we conclude 
in Sect.|5] 

2. NSE, thermal weak equilibrium and dynamic 
beta-equilibrium 

At high temperatures (T > 4 GK =s 0.3 MeV) the photon 
energy in an astrophysical plasma is high enough to dissoci- 
ate nuclei into nucleons. At the same time, for densities p w 
10 6 - 10 13 g/cm 3 , rapid nuclear reactions lead to the formation 
of nuclei. An equilibrium is reached between strong and electro- 
magnetic reactions. This is known as the NSE, where the abun- 
dances of nuclei follow simply from Saha equation for a given 
density, temperature, and electron fraction. Although this is well 
known, we briefly summarize the relevant equations. 

When NSE is reached, production and destruction of nuclei 
occur at the same rate and the chemical potentials of nuclei and 
nucleons have to satisfy: 

yU (Z ,A) = ZfJ. p + (A - Z)yU„ , (1) 

where the chemical potentials include the total mass 

Hi = niiC 2 + kTrji + \i c . (2) 



Here 77, is the degeneracy parameter, k the Boltzmann constant, T 
the temperature, and fx c the Coulomb contribution to the chemi- 
cal potential (see appendix in lJuodagalvis et al.. 2009). The elec- 
tron chemical potential is given by fi e - m e c 2 + rj e . 

The abundance T,- = n,/n of isotope i is given by the ratio 
of its number density n, over the total baryon number density n. 
The number density n,- is related to the degeneracy parameter r\\ 
assuming Maxwell-Boltzmann statistics by: 



A 3 



(3) 



with the partition function G,{T) = £(27* + l)e~ Ek/kT and the 
thermal wavelength A, — ^2nh 2 1 nijk T . For the partition func- 
tion we use values calculated bv lRauscherl d2~003). 

Protons also follow Maxwell-Boltzmann distributions over 
the densities of interest, but neutrons become degenerate for p > 
10 11 g/cm 3 (the transition from the non-degenerate to degenerate 
regime of course depends on the temperature). Therefore, we 
use Fermi-Dirac distributions for nucleons (both neutrons and 
protons), 

«i = ^^m^^m/zim^i) +#5 /2 fa,A)] , (4) 

where B = kT/mc 2 is the so-called relativistic parameter, which 
is almost zero for the range of temperatures considered here, and 
Tk are Fermi functions: 
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(5) 



The electron mole number is defined as Y e = n e /n, with the 
lepton number density n e = n e - - n e +. The electron, n e -, and 
positron, n e +, number densities are related to the respective 
chemical potentials by Eq. (0]l with the additional condition 

He = fie- = 'Me*- 

Combining Eqs. (Q} and (O, the abundance of the nucleus 
with Z, A is 



Y(Z,A) = 



J(Z,A) 



/ kTAm u 

p/m u \ lirh 2 



3/2 



(A-Z))7„ +Zr/ p B/kT Zfi CJ , -ii cS za> 



(6) 



Here p is the baryon density and m u the atomic mass unit. Note 
that small corrections (l ess than 1 %) arise from using m, = Am u 
dSeitenzahl et al.L 120091) . To calculate the NSE composition, we 
take the binding energi es, B - Zm p c 2 + (A - Z)m n c 2 — m^j^c 2 
from Aud i et al. j|20 03) and supplemented by theoretical values 
from Moller et al. (1995). In addition, the degeneracy parame- 
ters rj„, T] p are determined by charge neutrality and mass conser- 
vation, 



YjZiYi = Y e , 

= 1. 



(7) 
(8) 



Eqs. © - (O can be solved to obtain tj„, j] p , and Y(Z, A) for given 
temperature, density and electron fraction. 

The electron fraction is generally determined by weak pro- 
cess that are much slower than the strong and electromagnetic 
nuclear reactions responsible for maintaining NSE. However, 
if the system remains at a given temperature and density long 
enough, weak reactions can also reach equilibrium. If the neu- 
trinos produced by the weak interaction do not leave the system, 
the equilibrium reached is denoted as thermal weak equilibrium. 
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In order to determine the evolution of the system before thermal 
weak equilibri um is achieved neutrino transport calculations are 
necessary (see Murphv, 1980). Once, thermal weak equilibrium 
is achieved, the neutrinos have a Fermi-Dirac spectrum with a 
temperature, T, and neutrino chemical potential p v = —fiy. The 
state of the system is completely determined by the tempera- 
ture, density and lepton fraction (relative number of electrons 
and neutrinos). 

Once thermal weak equilibrium is reached, the composi- 
tion does not change unless the temperature or density varies. 
Therefore, entropy is conserve d which implies that TdS = 
2, yU/dF,- = (see MeygE I1993L for useful discussion of equi- 
librium in nucleosynthesis), thus 

flyYy + fi e Y e +^fliYi 
i 

= (Jl e - p v )Y e + (ji p ~H„)Yj Z ^ + V* Yj = ° ' (9) 

i i 

here the right-hand side is obtained from the assumption of NSE 
and the condition of lepton number conservation that implies, 
Y v - ~Y e , The time variation of the abundance of nucleus i is 
given by F;. Using Eqs. and © this reduces to 

Qi e + fi p - p. v - p. n )Y e = 0. (10) 

The value of Y e for which equilibrium is attained can be obtained 
from the thermal weak equilibrium condition 

p e +p p =li v +p n . (11) 

This implies that the value of the neutrino chemical potential is 
a function of lepton fraction, temperature, and baryon density. 
Alternatively, the equilibrium Y e can be obtained from 

Y e = - 2^(/le C + + 4 e + A 'inl3 ) Y i 
i 

+ £(4c + 4- + < + = 0, (12) 

i 

which requires the knowledge of the weak interaction rates for 
electron capture (ec), positron emission antineutrino ab- 
sorption (v e ), inverse of electron emission (in/? - ), positron cap- 
ture (pc), electron emission (J3~), neutrino absorption (v e ), and 
inverse of possitron emission (in/3 + ) for a large ensemble of nu- 
clei. Eqs. (fTTT i and (fT2l represent equivalent ways of obtaining 
the steady state value of Y e , if entropy is constant. Eq. (fTTT) is 
certainly advantageous as it depends only on bulk nuclear prop- 
erties like binding energies and partition functions. However, 
its applicability is rather limited because at stellar densities be- 
low ~ 10 11 g crrT 3 neutrinos scape and consequently a neutrino 
chemical potential cannot be defined. In these cases, a dynami- 
cal beta-equilibrium can still be achieved and its equilibrium Y e 
value is given by Eq. ( fT2l suppressing the neutrino absorption 
rates: 

Y e = - 2(4 + 4+) y <' + Z ( 4c + 4-) y ' = °- ( 13 > 

i i 

Nevertheless, we will show in the following that Eq. ( TTTb to- 
gether with the assumption p v = provides a very good approxi- 
mation to the solution of Eq. ( fT3l ), i.e., thermal weak equilibrium 
with fiy = reduces to dynamical beta-equilibrium. In this case 
the neutrino density depends only on temperature: 

_ 1 r p 2 

" V "" P ~ 2^3 J () l +e xp( P c/kT) P 

= 7.61 x 10 27 (— ) cm" 3 , (14) 



with p being the neutrino momentum and pc its energy. At 
low temperatures the neutrino densities are so small that can 
be neglected and consequently the condition fi Y = reduces to 
n v = 0. This is the situation in c old neutron sta rs where Eq. (fTTT i 
with p v — is commonly used dWebeit 1 1999b to determine the 
composition. However, at higher temperatures the assumption of 
/i v = results in neutrino densities comparable to electron den- 
sities and consequently thermal weak equilibrium, Eq. (fTTT i. and 
dynamical beta equilibrium, Eq. (fT3l . predict different values of 
Y e . Under these conditions, o ne has to solve Eq. (fT3l using a full 
set of weak interaction rates (lOdrz vwoleki 120091) . 

As discussed in the introduction, the determination of weak- 
interaction rates on nuclei is a very challenging nuclear struc- 
ture problem that requires sophisticated many body calculations. 
Reliable weak-interaction rates are necessary to determine the 
evolution of the system and t o account for the energy loss by 
neutrinos (Odrzywolek, 2009). However, in this paper we are 
not interested in the evolution towards equilibrium, but only in 
the value of the equilibrium electron fraction for a given tem- 
perature and density. Therefore, only a simple description of the 
weak interaction rates is necessary. 

We start considering the case of thermal weak equilibrium 
under the assumption that p v = 0. This has the advantage that we 
can use Eqs. (fTTT i and (fT2l to check the validity of our implemen- 
tation of the weak interaction rates. Once these are validated, we 
can can easily generalize our approach to the case with neutrinos 
leaving the system, by setting the neutrino densities to zero. Let 
us consider a couple of nuclei (A, Z) and (A, Z - 1) connected by 
the weak interactions: 



e~ + (A,Z) ?± (A,Z- l) + v e , (15a) 



(A,Z)<±(A,Z-l) + e + + v e , (15b) 



(A,Z-l)*±(A,Z) + e-+v e , (15c) 



e + + (A,Z-l)^±(A,Z) + v e . (15d) 

If neutrinos are present, these reactions can operate in both di- 
rections. Only the right direction is relevant for cases where neu- 
trinos leave the system. One can approximate these reactions as- 
suming that they occur by a transition from the ground state of 
the initial nucleus to the ground state of the final nucleus. This 
is an exact result for neutrons and protons and has been shown 
to provide a reasonable description of e lectron capture rates fo r 
a broad range of nuclei (see figure 2 of lLanganke et all 120031) . 
Here we generalize this prescription to the other rates in Eq. (TT3T >. 
Our aim is not to obtain accurate weak-interaction rates but an 
approach that reproduces the balance between forward and re- 
verse reactions once equilibrium is achieved. In this sense our re- 
sults should not be considered as a substitute to the different tab- 
ulations of weak interaction rates available in the literature. We 
use the following expressions to calculate the weak-interaction 
rates in Eq. (TT3T >: 
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= C^E\E-Qf 


fe(E,H e ) 
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{\-f p (E,-He)) 
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V 


= C^E\Q-Ef 


(l-f e (E, Me )) 




Anr 


= C C E 2 (E + Of 


f p (E, -n e ) 

(l-ME + Q, -fx v ))dE, 


(16d) 


K 


=cs:e\e- q ? 

J VV ec 


(l-f e (E,fi e )) 
f v (E-Q,fi v )dE, 


(16e) 


io/3* 


= c£ Q cl EH-Q-Ef 


f p (E, -n e )f v {-Q - E^ v )dE,{m 




= c C> eH q- e)2 


fe(E,Li e )MQ - E, -n v )dE, 


(16g) 




= cQe\E + Q) 2 


(l-f p (E,- Me )) 
ME + Q, -fx Y )dE, 


(16h) 



where the first four equations correspond to the forward direc- 
tion in Eq. $15[ describing the rates for electron capture (ec), 
positron emission (B + ), electron emission {B), and positron cap- 
ture (pc), respectively. The last four equations describe neutrino 
absorption (v e ), the inverse of positron emission (in/3 + ), the in- 
verse of electron emission (in/3 + ), and antineutrino absorption 
(v e ). Q is the transition Q value, i.e. Q = M(A,Z - l)c 2 - 
M(A,Z)c 2 , that is positive for neutron-rich nuclei and protons. 
M(A, Z) is the nuclear mass. The quantity C is given by 



C 



B\n2 



(17) 



K(m e c 2 ) 5 ' 

with K = 6144 s dHardv & Towneitl2009h and B the nuclear ma- 
trix element. We use t he value 5=1+ 3g 2 = 5.76 for nucleons 
and B = 4.6 for nuclei (Lang anke et al.Ll2003h . We have checked 
that our results are not sensitive to 20% variations of this value. 
/ is the Fermi-Dirac distribution 

f(E,fi) = — — — . (18) 

1 + exp[(£ - LQ/kT] 

The quantity w ec = max(m c c 2 , Q) is the threshold energy for 
electron capture, and w pc = max(m e c 2 , —Q) for positron capture 
rate. The integral in Eqs. (I16bl ) and ( I16fl ) are defined only for 
-Q > m e c 2 . Similarly, the integrals in Eqs. ( 116cl i and ( |16g| l are 
defined only for Q > m e c 2 . 

The value of the equilibrium electron fraction agrees up to 
3 significant figures when we solve Eq. ( fTTT i or Eq. < TT~2T >. This 
justifies our simple scheme, Eq. ( fTBI l, for calculating the weak 
interaction rates. 

Dynamic beta-equilibrium, Eq. ( fT3l . can be obtained from 
the rates in Eqs. ( fT6b setting the neutrino distributions to zero. 
The equilibrium electron fraction in dynamic beta-equilibrium is 
typically smaller than the one calculated in thermal weak equi- 
librium. As in the latter, neutrino absorption results in an in- 
crease of the electron fraction. 

3. Results 

We solve the NSE equations (Eqs. (|6]l-(|8]l) under the constraint 
of thermal weak equilibrium, Eq. ( fTTT i. assuming /j. v = for tem- 
peratures from 5 to 150 GK (0.43 - 13 MeV) and densities in the 
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Fig. 1. The electron fraction for NSE with thermal weak equi- 
librium is shown by color contours for a range of densities and 
temperatures where NSE can be applied and which are relevant 
in different astrophysical scenarios. The solid black line shows 
the core-collapse trajectory for a 15 M Q (A. Marek priv. com.); 
the dashed lines represent the presupernova evolution of a 15 M Q 
(black) and 25 M Q (blue) stars; the red dashed lines represent a 
type la trajectory (F. Ropke priv. com.); and the solid yellow 
lines correspond to the evolution of the proto-neutron star sur- 
face. Not all these scenarios will reach beta equilibrium and their 
electron fraction can therefore be different to the one shown here. 



range 10 6 - 10 12 g/cm 3 . This results in the equilibrium electron 
fraction shown in Fig.[T]and in the composition shown in Figs. [4] 
[5] and [6] Because only small changes are expected at high tem- 
peratures, we plot all quantities only over a smaller temperature 
range. The whole data table is available on request. 

In Fig.[T]the color contours show the electron fraction assum- 
ing thermal weak equilibrium and the lines represent trajecto- 
ries (density and temperature evolution) of various astrophysical 
scenarios: presupernova, core collapse, proto-neutron star sur- 
face, and type la supernovae. The collapse trajectory (A. Marek 
priv. c om.) of a 15 M sta r based on the presupernova mod- 
els of IWooslev et al.l (l2002f) is shown by the solid black line. 
When the inner core of a massive star collapses the evolution 
is much faster than weak reactions, therefore beta equilibrium 
is not achieved. Once densities above 10 11 g cirT 3 are reached 
neutrinos begin to be trapped and the assumption of fi v = is 
not valid any more. Therefore, the equilibrium Y e differs from 
the value obtained in detailed hydrodynamical simulations with 
neutrino transport. However, we include this trajectory in our 
figure for completeness, to have an idea of the range of temper- 
atures and densities we are s tudying. Presuperno va trajectories 
for 15 M Q and 25 M stars (IWooslev et al.Ll2002l) are shown by 
the dashed lines. This is discussed in detail in Sect. 14.11 The 
solid yellow line in the high temperature region represents the 
evolution of the hot proto-neutron star surface, where both ther- 
mal weak equilibrium and NSE are fulfilled (Sect. l4~2i i. And the 
red dashed line represents a trajectory of a type la supernova (F. 
Ropke priv. com.), where the evolution normally proceeds much 
faster than the weak interaction timescale and consequently no 
equilibrium is achieved (Sect. [43). 

For completeness we show in Fig.[2]the electron fraction ob- 
tained in NSE assuming dynamic beta-equilibrium, Eq. ( fT3l . for 
the same temperatures and densities as shown in Fig.Q] As dis- 
cussed in Sect. [2] we expect that dynamic beta-equilibrium and 
thermal weak equilibrium with jj. v = give similar electron frac- 
tion whenever the absorption of thermal neutrinos is negligible. 
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Fig. 2. The electron fraction for NSE with dynamic beta- 
equilibrium is shown by color contours for the same range as 
in Fig.Q] The black line represent the region where Ay e = Ap-. 



This occurs on the left side of the solid line shown in Fig. [2] This 
line marks the conditions for which the antineutrino absorption 
would become equal to beta decay if neutrinos did not escape. At 
high temperatures the results for both equilibria are not so differ- 
ent because the composition is dominated by neutrons and pro- 
tons and moreover positron captures also become important. In 
fact, at low densities (below the curve marked as r\ e — in Fig. [3} 
one has that p. e * p e + * implying that the contributions of 
electrons/positrons and antineutrinos/neutrinos are similar. This 
explains why both equilibria lead to same electron fraction for 
Big Bang conditions. The biggest differences between both ap- 
proaches appear in the region where the composition changes 
from heavy nuclei to nucleons (Fig. [3]). Moreover, it is in this 
region where the electron fraction drops more abruptly and con- 
sequently its value becomes rather sensitive to the equilibrium 
approach used. 

Although solution of the full equations gives the value of the 
equilibrium electron fraction, some simple approximations help 
to understand the results. We divide the p - T plane in differ- 
ent temperature and density ranges to explain the behaviour of 
the electron fraction. These ranges are shown in Fig. [3] The two 
horizontal lines correspond to electron and neutron degeneracy 
equal zero. We consider densities below the line of r\ e = to be 
low and above the line of r\ n = to be high. The temperatures 
are considered to be low when heavy nuclei or alpha particles are 
present, which in the figure correspond to the regions labeled as 
Xheavy > 0.5 and X a > 0.5, respectively. 

In NSE at high temperature (T > 15 GK) one can assume 
that only nucleons are present because all nuclei are dissociated, 
thus Y e = Y p and Y„ = 1 - Y e , For low densities conditions are 
similar to the Big Bang, where the composition consists mainly 
of non-degenerate neutrons, protons, and electrons. Therefore, 
we have p e — and p„ w p p that together with Eq. (O give: 



1 



1 + e-A/K" 



(19) 



where A = m„ - m p = 1 .293 Me V. Under these conditions Y e can 
reach rather high values when the temperature becomes smaller 
than A. At higher densities (p > 10 8 - 10 9 g/cm 3 ), electrons 
become degenerate and this approximation is not valid anymore. 
This happens when the density is above: 

p> 1.132- 10 1 — (0.6773 + 2.256 kT), (20) 




20 25 30 35 
Temperature [GK] 



50 



Fig. 3. Lines limit the regions where different approximations for 
Y e can be used. Horizontal lines indicate where the degeneracy 
of electrons (77^) and neutrons (77,,) is zero. The regions where 
nuclei and alpha particles dominate the composition are at low 
temperatures. 



ng rj e ss and using the val- 
73) 



which co rresponds to Eq. dUi assu ming 

ues from Blu dman & van Riperl d 19771) for the Fermi functions 
CF i/2 (t7 = 0) = 0.6773 and ^3/2(77 = 0) = 1.153). This line is 
shown in Fig.[3]labeled as rj e = 0. 

For densities above the line given by Eq. ( f20T > the electron 
fraction can be estimated by taking into account that the electron 
number density is approximately given by: 



(kT? 1 , 



(21) 



From this equation we get directly T] e . If we assume that the nu- 
cleons follow Maxwell-Boltzmann distributions (Eq. ©), their 
chemical potentials are: 



Pn(p) 



Vp) c 



kT\n 



A 



n{p) 



(22) 



Inserting all this in the beta equilibrium equation (p„ — p p + p e ) 
we obtain a non-linear equation for the electron fraction: 



kT\n 



3/2 , 
P I 1 



3i>\' /3 he _ Q 
87rm„ / kT 



(23) 



This can be solved numerically using as guess value the Y e of 
Eq. ( fT9l . At high density, low Y e (more neutron-rich material) is 
expected from Eq. (l23l . The last term in Eq. (l23l increases as the 
density increases, unless Y e decreases. Moreover, a reduction of 
Y e leads to an increase of the third term, which has opposite sign 
to the term containing the density. The exact solution, shown in 
Fig. Q] demonstrates the neutron richness of the matter at high 
densities. 

For low temperatures (T < 15 GK) the situation is more 
complicated because of the presence of bound nuclei. We esti- 
mate the densities and temperatures at which nuclei dominate th e 
composition following the derivation of lHovle & FowleJ (fl960). 
This leads to nuclei (represented by 56 Fe) corresponding to half 
of the mass at: 



logp = 11.62+ 1.5 log T 9 



39.17 

~tV' 



(24) 
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Fig. 4. Mass fraction of alpha particles for the same temperature 
and density range as in Fig. Q] In the white regions the mass 
fraction of alpha particles is negligible. 

In a similar way, taking mass fraction (X, = A,T,) of alpha par- 
ticles X a = 0.5, alpha particles represent half of the mass, for 
which one has: 

logp = 29.68 +4.5 log T 9 + log - , (25) 

where X a , X n , and X p are the alpha, neutron, and proton mass 
fractions, respectively, which are obtained from mass and charge 
conservation: 

1 = X a + X p + X„ , (26) 

Y e = -X a +X p . (27) 

At low densities, between the two lines given by Eqs. (1241 and 
d25l l. there are mainly alpha particles, therefore Y e - 0.5. Above 
the line given by Eq. ( l24l i. the composition is more complicated 
due to the presence of different nuclei. For high density the nu- 
clei become neutron rich. We can delimit the region where heavy 
nuclei are present using simple approximations, as described be- 
fore. However, the detail composition and thus the electron frac- 
tion in this region can only be obtained by numerically solving 
NSE equations under the assumption of thermal weak equilib- 
rium or dynamic beta-equilibrium, as described in Sect. [2] 

The mass fraction of alpha-particles is shown in Fig. [4] 
Alpha-particles appear mostly at low density and temperature, 
between the region where nucleons dominate the composition 
(high temperature) and the region where heavier nuclei are 
formed from the alpha particles. Notice that the formation alpha 
particle leads to a change in the behaviour of Y e at low densi- 
ties (Fig. [TJ. The mass fractions of heavy nuclei (A > 4) are 
shown in the panel (a) of Fig. [5] they dominate the composition 
at low temperatures and high densities. Figure[5]shows the aver- 
age mass, proton and neutron numbers of the nuclei in panels (b), 
(c), and (d) respectively. In addition we find that the abundances 
of light nuclei (A < 4) are significant at high temperatures and 
relative high densities, as shown in Fig. [6] Such conditions are 
found in the outer layer of the proto-neutron star which is still 
hot as shown by the trajectory (yellow line) in Fig.Q] 

4. Discussion and astrophysical implications 

There are several astrophysical environments where our results 
are interesting (Fig.[TJ. These are situations where sufficient time 
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Temperature [GK] 
Fig. 6. Mass fraction of light nuclei (A<4). 



elapses at high temperature for NSE to prevail and for weak in- 
teractions to approach equilibrium. Even when equilibrium is 
not fully achieved, our results give an interesting lower bound 
on Y e in cases where electron capture is dominant. 

4. 1 . Presupernova evolution 

We investigate here two presupernova models that were al- 
ready studied in de tail in the frame of weak interactions by 
Hege r et all (I2001 alibi) . Fig. [7] shows the results for the 15 M Q 
star in the left column and for the 25 M star in the right 
one. All results are plotted versus time before collapse which 
is at time zero. The dotted l ines are the Y e of the old models 
of IWooslev & Weaver! d 1995b (WW). These mod els were com- 
puted with electron-cap ture rate s of [Fuller et all d 19801) (FFN ) 
and beta-decay rates of iHanserJ dl968l) : iMazurek et al.l ([1974). 
The das hed lines were obtained when using th e weak-interaction 
rates of[Langanke & Ma rtmez-Pine do (2001) (LMP) in the cal- 
culation of the pres upernova models. The se two results were 
already presented in iHeger et all d2001albl) . where they discuss 
that differences in the electron fraction are mainly due to discrep- 
ancies between the old beta decay rates and the new ones. The 
bottom panels in Fig.[7]s how the electron-capture and beta-d ecay 
rates as calculated by Langanke & Martme z-Pinedo J (|2001l) and 
used fo r the presupernova models presented in iHeger et al.l 
d2001allbh . 

Before silicon burning, the evolution of the electron fraction 
is very similar for WW and LMP because electron captures dom- 
inate over beta decays (see bottom panels in Fig. |7]) and these 
rates are similar for FFN and LMP. In this first phase when Y e 
starts to decrease from the starting value of 0.5, the temperature 
is too low and the assumption of NSE is not valid and there are 
no values for the equilibrium electron fraction. 

During silicon shell burning the electron fraction continues 
dropping and beta decay becomes important and comparable 
to electron capture. For the 15 M star b oth rates are equal 
and a dynamic beta-equil ibrium is reached ( Aufde rheide et all 
Il994at iHeeer et aU 1200 lbl) . After silicon burning the electron 
fraction stays almost constant although electron captures be- 
come dominant over beta decays because the dynamical evolu- 
tion is faster than the cap ture rates. Therefore, as it was pointed 
out bv lHeger et al.l d2001bl) . the final electron fraction is set be- 
fore the core begins its final contraction. The most important 
period for determining the core structure is thus the silicon shell 
burning phase. 
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Fig. 5. The panel (a) shows the mass fraction of heavy nuclei (A > 4), in the white region no heavy nuclei are present. The other 
panels give the average mass (b), atomic (c), and neutron (d) numbers. The blue line in the mass fraction plot marks the contour 
where heavy nuclei are half of the mass. In the average number plots, the black lines correspond to magic numbers and the red ones 
to 56 Fe. 



We use the density and temper ature evolutions shown in 
Fig. [T]for the presupernova models of Woosl ev"& Weaver! (l 995 ) 
and calculate the equilibrium electron fraction. This Y e corre- 
spond to the solid lines in the upper panels of Fig. [7] Notice that 
there are two solid lines, the black one results from the assump- 
tion of thermal weak equilibrium, Eq. ( fTTT i. while the green one 
is obtained under the assumption of dynamic beta-equilibrium, 
Eq. (fT~3T >. The resulting Y e values are almost identical. For the 
15 M Q star, where dynamic beta-equilibrium is achieved in the 
simulation (see lower right panel of Fig. [7J, all the Y e values are 
almost identical for times around 10 4 s. For the 25 M star, dy- 
namic beta-equilibrium is never reached and consequently our 
equilibrium Y e is slightly smaller. 

Our equilibrium electron fraction represents a static equi- 
librium and is obtained independently of the weak-interaction 
rates. Such an equilibrium is reached if the temperature and den- 
sity are kept constant long enough. There are two points that a 
set of weak-interaction rates has to fulfil to be used in presu- 
pernova models: 1) if dynamic beta-equilibrium is reached, the 
electron fraction should be very similar to our equilibrium elec- 
tron fraction; 2) the electron fraction in the presupernova models 
can not drop below the e quilibrium value we calcula te. The older 
presupernova models o f IWooslev & Weaverl d 19951) do not ful- 
fil these two constraints. Auf derheide et al.h 1994al) : lHeger et ail 
(2001b) suggested that older models needed to be recomputed 
with consistent set o f weak-interaction rates. The new results 
dHeger et all 1200 lb) based on better determination of the rates 
satisfy the two constraints. This test should also apply to future 



compilations of weak-interaction rates and their implementation 
in presupernova simulations. 

4.2. Proto-neutron star 

Iron core collapse marks the end of the presupernova phase. 
When the collapse starts the dynamical evolution is very fast 
compared to the weak interaction timescale, thus the assumption 
of equilibrium is not valid. But when collapse continues, neutri- 
nos are trapped due to the increasingly higher central densities, 
and thermal weak equilibrium can be reached again. However, 
neutrinos are trapped and their chemical potentials are not negli- 
gible anymore (fj. v + 0). Therefore, our results give just a rough 
estimation of the electron fraction during collapse, but neverthe- 
less are useful to have an idea of the composition in the collaps- 
ing material. In Fig.Q]the trajectory of the collapse is in a region 
where almost only heavy neutron-rich nuclei and neutrons are 
present (Fig. O. The amount of alpha particles (Fig. |4j and light 
nuclei (Fig. [6) is negligible. 

When central density gets above the density of nuclear mat- 
ter, collapse stops and falling matter bounces back leading to the 
formation of a shock. The way this shock propagates through 
the outer layers of the sta r in a supernova e xplosion is still a 
problem under discussion dJanka et all 120071) . But it is known 
that a compact object forms in the center and cools by neutrino 
emission. Neutrinos start to decouple from matter at the neu- 
trinosphere that is placed in the outer layers of the newly born 
proto-neutron star. In this region we assume still that [i v m 0. 
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Fig. 7. The electron fraction (upper plots), the electron-capture, and beta-decay rates (bottom plots) as a function of time till collapse 
are shown for two presupernova models: 15 M (left) and 25 M Q (right). In the electron fraction plots the dotted lines (blue) and 
dashed lines (red) were already shown b y iHeeer et al.l d2001al) and correspond to the old models of IWooslev & Weaverl d 19951) and 
to the new ones of He ger et alj d2001bl) . respectively. Here we have added the solid lines, which results from assuming NSE and 
thermal weak equilibrium (solid black line) or dynamic beta-equilibrium (solid green line). 



Also NSE is satisfied because the temperatures are high enough, 
T ~ 50 ■ 10 9 K. We compare our equilibrium electron fraction 
with the one obtaine d in hydrodynamical simulations of core- 
collapse supernovae dArcones et all 120071) . Figure [8] shows the 
post-bounce evolution of the electron fraction. T he dashed line 
is the result of the hydrodynamical simulation of lArcones et alj 
(2007) for their standard case (M1511rl), and solid line is the 
equilibrium electron fraction computed for same density and 
temperature evolution as for the dashed line. There is a clear dis- 
crepancy between the two electron fractions that increases with 
time. This disagreement comes from t he dif ference in the com- 
position (see Fig. 4 in lArcones et all 120081) . The equations of 
state that are generally used in supernova simulations include 
only nucleons, alpha particles, and a representative heavy nu- 
cleus. In our calculation of NSE with beta equilibrium, we take 
into account all nuclei. Figure [6] shows that the abundances of 
light elements (A < 4) are not negligible in the outer layers 
of the proto-neutro n star. Similar results are f ound with other 
EoS, for example [Horowitz & Schwenkl d2006l) and lTvpel et alj 
(2010). Light nuclei are formed at expenses of the few free pro- 
tons. This produces also a reduction of the antineutrino absorp- 
tion and electron captures, which leads to an increase of the elec- 
tron fraction compared to the case where only neutrons and pro- 
tons are present (see Eq. dT2b). 

These light nuclei have an impact on the prope rties of neutri- 
nos em it ted during the first se c onds a fter bounce dArcones et alj 
d2008l) 1. ISumivoshi & Ropkd d2008l) have also pointed out the 



influence of the light elements during the post-bounce evolution. 
These calculations were both done in a post-processing step, 
therefore the real effect of light elements in the dynamics of the 
supernova evolution is still an open issue that has to be studied 
in detail by using an equation of state that includes light nuclei 
in supernova simulations. 

4.3. Type la supernovae and Accretion-induced collapse 

Other astrophysical scenarios where the electron fraction plays 
an important role involve the evolution of white dwarfs in bi- 
nary systems. A white dwarf accreting up to the Chandrasekhar 
limit can have two possible outcomes. If the density and tem- 
perature become sufficiently high to ignite explosive nuclear 
burning, then the white dwarf may become a type la super- 
nova. But it could also happen that electron capture in the 
ash reduces central temperature and pressure, and the white 
dwarf instead collapses to a neutron star. This lat ter outcome 
is known as "accretion-induced collapse", or AIC dCanal et all 
1990; IWooslev & BaronL Il992t iFrveret all Il999h . Which out- 
come occurs depends on the ignition density which is in turn 
determined by the initial white dwarf properties (ma ss, compo- 
sition , and temperature) and on the accretion rate dCanal et all 
IT99Qh . 

The outcome also depends on the weak interaction rates em- 
ployed in the study and the multi-dimensional character of the 
explosion. Burning raises the temperature at constant pressure 
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Fig. 8. The evolution of the electron fraction at the proto-neutron 
star surface is shown. The solid line is obtained by assuming 
NSE and thermal weak equilibrium, while the dashed line results 
from the hydrodynamical simulation of the sup ernova explo- 
sion a nd the subsequent neutrino-driven wind of lArcones et al.l 
(2007). 



and, if Y e is constant, that must reduce the density. The density 
inversion behind the flame drives a Rayleigh-Taylor instability 
that greatly enhances the advance of the burning front. However, 
since the pressure is mostly due to relativistic electrons, a de- 
crease in Y e due to electron capture in the ash means that the 
post-burning density will be larger (since pressure depends on 
the product pY e ). If Y e decreases sufficiently, the density beneath 
the flame may even rise, thus suppressing the Rayle igh-Taylor 
instability (see Sect. 4 in Timme s & Wooslevl 119921) . If Y e de- 
creases in a sufficient fraction of the white dwarf before the (now 
slowly moving) flame makes its way out to regions where the 
density is lower and electron capture is reduced, the white dwarf 
collapses. At densities approaching 10 10 g ctrr 3 electron capture 
is so efficient that it would drive Y e to very low values were it not 
for the inhibiting effect of beta-decay. Yet there are, so far, no ac- 
curate calculations of beta-decay rates for Y e below about 0.41. 
Having a reliable lower bound to Y e at a given density will thus 
be useful in future simulations. 

During this phase of the white dwarf evolution, tempera- 
tures are high and densities relatively low. Consequently, the 
produced neutrinos escape. Under such conditions, the lower 
limit of the electron fraction is thus obtained from dynamic 
beta-equilibrium, Eq. ( fT3l ). If thermal weak equilibrium were as- 
sumed, electron captures would compete with antineutrino ab- 
sorptions instead of p~ decay, as expected when neutrinos es- 
cape. Figure [9] shows the dynamic beta-equilibrium electron 
fraction for the range of temperatures and densities relevant for 
AIC and type la supernovae. In Tab. Q] some numerical values of 
the Y e in the figure are given (a complete table is available on 
request). 

The behaviour of the electron fraction is very different at 
high and low temperatures (Fig. [9]). At low temperatures the 
electron fraction decreases continuously with density. While a 
high temperatures (T > 12 GK) first decreases and then in- 
creases. This seems unintuitive but it can be understood con- 
sidering the change in composition. At low densities the compo- 
sition is neutrons, protons, and alpha-particles. With increasing 
density heavy nuclei are present and their fT decay produces an 
increase of Y e . At even higher densities, electron capture domi- 
nates and Y e decreases again. 



Ye 
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Fig. 9. The dynamic beta-equilibrium electron fraction (Eq. ( fT3l l) 
is shown for the range of temperatures and densities relevant in 
type la supernovae and AICs. In Table[T]numerical values of the 
equilibrium electron fraction are given. 

Finally, we notice that the electron fraction and thus weak- 
interaction rates are important for determining the composition 
of the ejecta in typical type la supernova. In these events the elec- 
tron captures are very important during the flame propagation 
because they lead to a reduction of the electron fraction which 
is the parameter that controls the isotopic composition of the 
ejecta. 

5. Conclusions 

Weak-interaction rates determine the evolution of the electron 
fraction, which is a key parameter affecting the composition and 
dynamics of stars in the late stages of stellar evolution and su- 
pernovae. Here we have presented an interesting limiting case 
that can be used to test the implementation of theoretical weak- 
interaction rates in stellar evolution simulations and to place 
lower bounds on the Y e in the explosion of white dwarfs. Our 
results depend only on chemical potentials and therefore are in- 
dependent of the weak-interaction rates. 

These equilibrium electron fractions can be compared to 
those found in presupernova models, where we ak-interaction 
rates play an important role dHeger et all l2001al) . The presu- 
pernova phase precedes the collapse of massive stars, therefore 
temperatures get high enough that NSE is reached, but the dy- 
namical evolution is too fast in most cases for beta equilibrium to 
be achieved. We find a lower limit for electr on fraction and show 
that th e Y e in the old presupernova models o flWooslev & Weaver! 
(119951) drops be l ow it. I n these presupernova models beta decay 
rates of lHansenl (Il968l) : iMazurek et all (11974 were used. These 
rates are too low compared to FFN and LMR Similar model s cal- 
culated using more recent weak rates (Heg eret"aTl l2001b) lead 
to an electron fraction that stays always above or equal to our 
lower limit. Moreover, for the 15 M G model where the rates pre- 
dict dynamic beta-equilibrium, we find that their electron frac- 
tion agrees with the equilibrium value of our calculations. 

The accretion onto white dwarfs approaching the 
Chandrasekhar mass could lead to type la supernovae or 
to AIC. The final outcome depends on the c ompetition betwee n 
explosive ignition and electron captures (ICanal et all Q990). 
Therefore, the electron fraction is a key parameter. In this envi- 
ronment densities and temperatures drive the NSE composition 
towards neutron-rich nuclei whose weak-interaction rates are 
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0.388 
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T 9 


= 6.010 


0.438 
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Table 1. Dynamic beta-equilibrium electron fraction for density values given in the first row and temperature values in the first 
column. The temperature is in units of 10 9 K and the density in 10 9 g cirT 3 . These temperature and density intervals correspond to 
theFig.|9] 



not yet known from theoretical models. Here we give a reliable 
lower limit for Y e , that can be used in simulations (data are 
available on request). 

In addition, our calculations show that the amount of light 
elements (A < 4) in the outer layer of the proto-neutron star is 
not negligible. After supernova explosion a hot proto-neutron 
star is born and cools emitting neutrinos. The surface of the 
proto-neutron star, where neutrinos decouple from matter, is hot 
enough to assume NSE and its evolution is rather slow, thus 
beta equilibrium is also satisfied. Changes in the composition of 
this region have an imp act on neutrino prop erties that could af- 
fect the nucleosynthesis (lArcones et afl.l2008l) and the explosion 
dSumivoshi & Ropkell2008l) . 
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